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l/^ ' Abstract 
(N ■ 

We extend the analysis of the Backgammon model to an ensemble 
with a fixed number of balls and a fluctuating number of boxes. In this 
ensemble the model exhibits a first order phase transition analogous to 

■ the one in higher dimensional simplicial gravity. The transition relies 
i on a kinematic condensation and reflects a crisis of the integration 
I measure which is probably a part of the more general problem with 

■ the measure for functional integration over higher (d > 2) dimensional 
5^ , Riemannian structures. 
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Recently we have proposed a scenario for the phase transition in higher 
Oh. dimensional simplicial quantum gravity based on the balls-in-boxes model 

r-| ! 0] also called Backgammon model in a different context It has been 

established that simplicial gravity has two phases: crumpled and elongated. 
In the elongated phase it is effectively described by the Branched Polymer 
model IQ. In the crumpled phase is characterized by the appearance of a 
singular vertices which gather around them an extensive part of the volume 
|]. The transition between those two phases is of first order 0, §]. The 
constrained-mean-field interpretation based on Backgammon model captures 
all those features except the order of the transition. 

In this letter we show that the modified balls-in-boxes model undergoes 
a first order phase transition of the same type as in simplicial gravity. 

We start with a brief review of the properties the Backgammon model. 
The model describes N balls distributed in M boxes. The partition function 
reads §, i: 

Zm,N = J2 Pi(ll)---Pi(lM)^qi + -+qM,N (1) 

The numbers of balls qi and qj in two boxes i and j are independent as 
reflected by the factorization of the total thermodynamic weight to one-box 



1 



weights p{q)^s. The factorization is weakly broken by the constrain on the 
total number of balls. This constrain introduces correlations between boxes 
and may result in the condensation of balls into one box 0]. 

We are interested in the limit of — *• oo and r = M/N fixed. We call 
the quantity r curvature. In this limit one expects the partition function to 
behave a^ 

ZM,^~exp(iV/(r) + ---) (2) 

where the dots stand for corrections growing slower than linearly in N. The 
function /(r) is the free energy density per ball. The quantity r is bounded 
from below by since the system has to have at least one box. Addition- 
ally to fix attention we introduce an upper bound on r by requiring that each 
box must contain at least one ball. These give together : < r < 1. 

The phase structure of the model is encoded in the analytic properties 
of the function (the log does not change those properties and is added for 
convenience) : 

K{fi) = -\ogJ2piq)e~'' (3) 
q 

In particular if the series under the logarithm has a finite radius of conver- 
gence then the model given by the partition function (|T]) has a very interesting 
phase structure P, |]. Denote the real value of at the radius of conver- 
gence by ficr, and restrict ourselves to the situation when the value of the 
derivative K'{ficr) is finite as happens for instance for the weights of the form 
p{(l) = g~"e~^'"''' when a > 2. In this case it can be shown that the system 
has two phases depending on the value of curvature r 0. 

The free energy density /(r) (0) can be found by the steepest descent 
method. One obtains the following expression : 

{-rK{fj,^{r)) + yU*(r) for r > r^r 
(4) 
— TKcr + /^cr for T < Tcr 

where yU*(r) is given by the solution of the saddle point equation 

- = ^'(^*) (5) 
r 

^One might alternatively formulate the limit in terms of M and p — N/M : Zm,n ~ 
exp (Mf(p) + • • •) in In fact this would look more natural, since p has an interpre- 
tation of the average number of balls per box and not the number of boxes per ball as 
defined by r. Moreover it would be more natural to keep M in front of / since then / had 
an interpretation of the free energy density per box. We use the convention with r and N 
on purpose since it is better suited for calculations in the ensemble with varying M which 
is the subject of the letter. 



2 



The positivity of the weights p{q)^s imphes that the function on the right 
hand side of (|^) is monotonic. From this follows that the solution /i^, is 
unique and that it is a monotonic real function of r. For r = 1 it is infinite. 
When r decreases so does /i^,. The equation has a solution for as long as r 
is larger then a critical value rcr- The critical value is obtained from the 
equation (||) when on the right hand side /x* = ficr reaches the singularity of 
K. At the critical value the free energy is given by the formula in the lower 
line of (^) where n^r = K{^cr) and sticks to this value even when one r is 
lower since /i* may not move beyond ficr- At this point the system undergoes 
a continuous phase transition g]. 

In simplicial gravity the analog of the quantity r is the ratio Nq/N^ where 
the Nq is the number of vertices and A^4 is the number of 4-simplices. The 
vertices play the role of boxes and the 4-simplices that of balls. The number 
of balls in one box is analogue of the number of vertices sharing one vertex 
ie vertex order. The sum of vertex orders Qi is constrained by the relation: 
qi + - ■ ■ + <iNo = 5N4. The ratio N0/N4 is a hnear function of Regge curvature. 
The average curvature is varied through a coupling to the number of vertices. 

In this letter we study the model with fluctuating number of boxes which 
more closely resembles the simplicial gravity scenario. The partition function 
of this model is : 

N 

Z{k,N)= J2 ZM,Ne^'' (6) 

M=l 

Substituting Zm,n by (|]) we get in the large N limit : 

Z{k,N) = N £drexpN(^Kr + f{r)) (7) 

where /(r) has the form (^). The integrand is the distribution of the cur- 
vature for a given k. For large N one expects the integrand to approach 
gaussian shape centered around being the solution of the saddle point 
equation : 

«: + /'(n)=0 (8) 

For K > Kcr this equation reduces to k = K{fisp{r^)) which has a unique 
solution for r^,. The value of r^, being the position of the gaussian distribution 
is the average curvature in the limit ^ 00. This situation continues as 

long as K > Kcr- 

Below the critical point k < Hcr the free energy /(r) is linear in r and 
the saddle point equation (H) has no solution. The integrand is not gaussian 
anymore but a decreasing function of r. It is exponential : expN{K — Kcr)r 
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Figure 1: Average curvature as the function of k, for infinite N (bold line) 32, 
64, 128, 256, 512 and 1024 balls 

for r < Tcr- For large only this exponential part matters in the integral 
(J^) and one finds that < r >~ 1/A^. 

Exactly at the critical point n = Kcr, the equation (H) is fulfilled by all r 
between and Vcr as can be seen from the second line of (^. The intuition 
behind this is that the curvature r must abruptly move between two separated 
regimes, namely gaussian at r = and exponential at r = 0. Therefore r 
stays undetermined by the saddle point equation. To fix the shape of the 
integrand one has to consider finite N systems. 

In the figure 1 we show average curvature (r) as a function of k, for weights 
p{q) = q^2. For k > it is a solution of the saddle point equation (||). It 
stops at and falls to zero. In the same figure we plot average curvature for 
some finite values of A^. The results were obtained by a recursive technique 
described in [Q. In the gaussian phase k > Kcr the curves lie close to each 
other indicating that the finite size effects are small there. For large k the 
curves approach asymptotically the upper kinematic bound r = 1. In the 
exponential (kinematic) phase k < Kcr finite size effects are stronger reflecting 
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the size dependence of the lower kinematic bound This also agrees 

with the results of numerical simulations of simplicial gravity, but there the 
kinematic bound is I/Vn rather than 

The steepest part of the curves corresponds to the pseudocritical region. 
It is a crossover region where the two phases coexist. The spread of the 
pseudocritcial region is 6k, ~ and reduces to one point for infinite N. 
In the figure 2 we show the distributions of r for = 512 and = 1024. 
One sees coexistence of two phases. The fraction of either phase depends on 
a particular choice of k in the pseudocritical region. We picked its values 
in such a way as to have more or less the same heights of both the peaks. 
The similar two peak structure has been observed in simplicial gravity for 
systems with 32000 and 64000 simplices 0, g 

This nicely extends the results of [Q showing that all the thermodynam- 
ical (non-geometrical) features of the simplicial gravity can be explained by 
this very simple model. It seems that the mean field sets in already for 
d > 3. The phase transition from Branched Polymer phase to the crumpled 
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phase is kinematical ie it is associated with system approaching the phase 
space boundary rather then diverging correlations (interactions). We think 
that this reflects the crisis of the measure in higher dimensional geometries. 
However it may also be that the gravity does not exist alone and needs some 
matter fields which would locally smooth the manifold and open a physical 
window. Unfortunately, fermions on random lattice seem to be for the time 
being far beyond numerical access. 

Of course the presented model can only capture the bulk features of the 
system as all the geometry has been integrated out, so its simplicity does 
not exclude the possibility of interesting geometrical behavior. One simple 
geometrical realization of Bacgammon is BP model [@]. In this model exactly 
at the phase transition system enters a third phase which has a different 



Hausdorff dimension p|, |T0[. It may be that the same phenomenon happens 
in simplicial gravity. 

The authors thanks Desmond Johnston and Jerzy Jurkiewicz for helpful 
comments and discussion. This work was partially supported by KBN grants 
2P03B19609. 
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